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Abstract. A search for the near- infrared water-ice absorption band was made in a number of very red OH/IR 
stars which are known to exhibit the 10/im silicate absorption. As a by-product, accurate positions of these highly 
reddened objects are obtained. We derived a dust mass loss rate for each object by modelling the spectral energy 
distribution and the gas mass loss rate by solving the equation of motion for the dust drag wind. The derived 
mass loss rates show a strong correlation with the silicate optical depth as well as that of the water-ice. The stars 
have a high mass loss rate (> 10 -4 Mq yr _1 ) with an average gas-to-dust mass ratio of 110. In objects which 
show the 3.1/tm water-ice absorption, the near-IR slope is much steeper than those with no water-ice. Comparison 
between our calculated mass loss rates and those derived from OH and CO observations indicates that these stars 
have recently increased their mass loss rates. 
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1. Introduction 

In order to estimate the dust mass loss rate in O-rich 
Asymptotic Giant Branch (AGB) stars, it is common to 
use the strength of the 10/im silicate emission. There 
seems to be a direct correlation between the 10/im excess 
emission and the derived mass loss rate estimated by CO 
(Skinner & Whitmore I1998f) . However, this does not ap- 

' ply when the silicate feature becomes self-absorbed or in 
absorption as mass loss rate increases. On the contrary, it 

' has been known that extreme OH/IR stars with strong 
silicate absorption show a very weak CO line strength 
which indicates low derived gas mass loss rate (Heske et 
al. ^990). Justtanont et al. 1)1996(1 investigated one such 
star, OH26.5+0.6 and found that the discrepancy can be 
explained by the fact that the star has recently undergone 
a superwind where the mass loss rate increases by two or- 
ders of magnitude. The low rotational transitions of CO 
mainly trace the cool gas mass loss rate in the outer part of 
the envelope, which has not been affected by this mass loss 
increase. The interferometric CO J=l-0 map using BIMA 
confirms the envelope size is marginally resolved (Fong et 
al. E302), consistent with a compact CO envelope due to 
the superwind. 

Accompanied by the deep 10/im silicate absorption, 
the near-infrared spectrum of some extreme OH/IR stars 
shows the absorption at 3.1/tm due to water-ice. Meyer et 
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al. I|1998fl suggested a correlation between the mass loss 
rate and the column density of water-ice in a few sample 
of stars observed. So far, there has not been a systematic 
search for the 3.1/xm absorption in these stars as a long in- 
tegration is needed in order to detect the absorption band 
when the infrared flux is already very low. A medium reso- 
lution spectrum is needed in order to discriminate between 
the type of ice condensed, i.e., amorphous or crystalline. 
Maldoni et al. (2003) suggested a significant fraction of 
water-ice in crystalline form condensed in the outflow of 
OH32.8-0.3. The crystalline ice has sharp features at 3.1, 
44 and 62/im, the latter two were first detected by the 
observation by the Kuiper Airborne Observatory (Omont 
et al. I1990fl . The libration band at 12/tm is more diffi- 
cult to discern due to the fact that it blends with the 
strong silicate absorption. Observations of OH/IR stars 
using the Infrared Space Observatory (ISO) Short- and 
Long Wavelength Spectrometers (SWS and LWS) show 
that these stars exhibit signatures of crystalline silicate 
dust (e.g., Cami et al. 1998), as well as water-ice (Sylvester 
et al. [TM% . 

2. Observations 

A sample of OH/IR stars were selected based on the deep 
10/tm silicate absorption, mainly from the IRAS low res- 
olution spectra (LRS) and the OH maser observations 
by Sevenster i|2002J) . Our observations were done using 
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Table 1. Observed fluxes of the OH/IR stars using SIRCA. 
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the Stockholm InfraRed CAmera (SIRCA) on the Nordic 
Optical Telescope (NOT) during July 22-23, 2003. We 
observed a total of 17 stars in the nodding and chop- 
ping mode. The stars were observed in the K, L', nar- 
row (AA = 0.15^m)- and wide-H 2 0-ice (AA = 0.3^im) 
bands. The central wavelength of the water-ice bands is 
at 3.07/im. The field-of-view is relatively large, 1'15" by 
1'15". Each star is flat fielded then flux calibrated using 
standard stars taken during the observation. A summary 
of the observations is given in Table Q 

Since we have images in several bands, it is a sim- 
ple matter to identify the OH/IR stars and determine 
their positions with the aid of the 2MASS images. The 
images at the K band in our observations are compared 
to the published map. In some cases, a star is not visi- 
ble in the K band so an interpolation of the positions of 
nearby stars was performed. Many of the fields show mul- 
tiple point sources and with the aid of the L' band, most 
can be distinguished as the OH/IR stars. For the field of 
star IRAS22177+5936, there is one very bright star which 
may be the target as the flux increases from K to L' bands. 
The measured flux in the narrow band is greater than the 
wide band at 3/xm which is most likely due to the satura- 
tion problem. We maintain this star in our sample as we 
have CGS4 and archive ISO data. From the difference of 
the wide- and narrow-water-ice band fluxes, we put the 
detection limit of the water-ice at 20% of the difference 
over the wide-band flux. We chose this limit based on the 
star OH26.5+0.6 which was shown to have water-ice ab- 
sorption at 3.1/im (Meyer et al. [1998). 

As a complement to our data set, earlier observations 
of six OH/IR stars done at the 3.8-metre UKIRT, using 
the Cooled Grating Spectrometer (CGS4) in August 1999 
are included to enlarge the SIRCA sample. Some of the 
stars overlap with the SIRCA sample (Table EJ). For most 



Table 2. A list of OH/IR stars using CGS4 and their 
observed positions. 
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stars, we have a spectral coverage between 2 to 4/im, punc- 
tuated by the earth's atmospheric absorption. The spectra 
were wavelength calibrated using arc spectra and then flux 
calibrated with standard stars. The absolute flux level is 
further adjusted to coincide with known photometry from 
either published data or SIRCA observations. This pro- 
cedure is done simply because the positions of these stars 
were not well known and the flux in the K band is very low. 
Some stars are mispointcd by several arcsec which means 
the flux level can be much lower than expected. This can 
be seen by comparing the coordinates of stars which are 
observed by both SIRCA and CGS4. Also, these stars are 
known to be variable with a period of up to 1 000 days. 

3. Determination of mass loss rates 

In order to derive the dust mass loss rate from our sample 
stars, a model fit to the spectral energy distribution (SED) 
of individual star is performed. This, in turn, is used to 
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Fig. 1. SED fit to the stars which do not show ice absorption. The symbols in the plots are x for the IRAS LRS, 
filled circles are the IRA S photometry, open circles are the SIRCA photometry and plus signs are CGS4 spectra. ISO 
spectra are shown as heavy solid line. 



calculate the total mass loss rate assuming the momentum 
driven wind. 



in the modelling as the major dust component is silicates 
which have absorption features in the lOyum window. 



3.1. Dust mass loss rates 

We derive the dust mass loss rate by fitting the SED for 
each star (Justtanont & Tielens ll992fl . We extracted IRAS 
LRS spectra which cover the wavelength range between 7 
to 24/im, along with the photometry from the point source 
catalogue. Where available, we also extracted archived 
ISO-SWS spectra. These complementary data are crucial 



The radiative transfer calculation used to model the 
stars is based on the work by Haisch (^979 ) . It assumes a 
spherically symmetric outflow, with a constant mass loss 
rate and outflow velocity. The thermal emission and scat- 
tering is calculated at a specified radial grid point and the 
emergent flux and the temperature profile of the dust cal- 
culated. We assume a grain size distribution for the dust 
to be the same as that of the interstellar medium (Mathis 
et al. I1977|l . As the distance to these stars are unknown, 
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Fig. 2. SED fit to the stars which show ice absorption. The 
same symbols are used for photometric data and IRAS 
LRS. The model with pure silicate component is shown 
as solid line while the composite silicate core-ice mantle is 
shown as dashed line. 

we assume a luminosity of 10 4 L Q and vary the distance 
until the overall spectrum fits. 

We use two different grain models to fit our SED in the 
envelope which shows water-ice absorption, namely, pure 
silicate grains and silicate core- water ice mantle. The dust 
parameters used to fit the absorption features are silicates 
based on OH26. 5+0.6 model by Justtanont et al. {HXty 
and crystalline water-ice from Schmitt et al. I|1998J) . We 
select the crystalline form of water-ice because there is ev- 
idence that the grains condensed in such an outflow show 
sharp structure at 3.1/im (Maldoni et al. I2003|) and on the 
bases of the far-IR features due to crystalline water-ice in 
the ISO spectra (e.g., Sylvester et al. ©99; Dijkstra et al. 
2003b ) The inner radius of dust shell is set to a point where 



the temperature is below the amorphous silicate conden- 
sation temperature of 1 000 K while the silicate core-ice 
mantle grains have a maximum temperature of 100 K, 
i.e., the temperature which water-ice condenses onto ex- 
isting silicate grains. The terminal velocity of the outflow 
of each star is well determined by either CO rotational 
line or via OH maser observations (see Table |3J). 

Figs, m [3 and show the model fits to individual ob- 
jects. Half of our sample of extreme OH/IR stars do not 
show convincing evidence for the water-ice absorption, de- 
spite the fact that the mass loss rate is very high. The 
derived mass loss rate and other parameters are listed in 
Table The range of the dust mass loss rate of these ex- 
treme OH/IR stars are between (1-4) 10~ 6 M Q yr -1 . In 
deriving the dust mass loss rate, we assume a constant 
dust velocity for all the grain sizes. Although this is the 
case for the smallest dust grains, the drift velocity for the 
largest size (0.25/zm in our model) can be up to 3 km 
s~ 4 . We address this issue of the dust velocity and how it 
affects the mass loss rate estimate in the next section. 



3.2. The dynamical mass loss rate 

In order to determine the total mass loss rate from a star, 
we used the derived dust mass loss rate and assume that 
the radiation pressure on the dust grains is the main driv- 
ing force of the circumstellar wind. The dust momentum 
is transferred to the gas and knowing the terminal velocity 
of the gas for each star, we can calculate the dynamical 
mass loss rate (Goldreich & Scoville fl97fi|l . The only un- 
known in solving the equation of motion is the mass of 
the star. We assume that the progenitor stars of these 
extreme OH/IR stars are more massive than the hot bot- 
tom burning limit of 4 Mq, otherwise they would end up 
as C-rich AGB stars. We assigned an arbitrary mass of 5 
Mq for all our stars. However, this parameter has a very 
small effect on the derived mass loss rate. A change from 
a 1 to 5 M© gravitational potential resulted in a decrease 
in the derived mass loss rate of < 30%. We also calcu- 
lated the dust velocity and fed it back to the calculation 
for the SED hence the value of the dust mass loss rate in 
Table 13 is consistent with the calculated dust velocity for 
the largest grain size. 

The main uncertainty in the derived dust and dynam- 
ical mass loss rates lie in the estimate of the distance to 
the star. In many cases, this is not known and the distance 
was estimated by assuming the stellar luminosity of 10 4 
L Q . It was shown that semiregular variables have a mean 
luminosity of 4 200L Q (Olofsson et al. 2002 ) and OH/IR 
stars are the more luminous counter part of these stars. 
An error in the distance estimate reflects directly the error 
in the mass loss rate. Fortunately, there is a direct scale 
between the mass loss rate and the distance, i.e., Md/D is 
constant in order to fit the same spectrum. Our derived 
total mass loss rates for these objects are very high, a few 
10~ 4 Mq yr" 1 , i.e., they are in a supcrwind phase. Such 
a high rate can be sustained by the high-mass end of the 
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Fig. 3. SED fit to the stars which show ice absorption. The 
lower than expected near-IR flux and needs a modification 
near-IR spectrum of OH 32.8+0.3 is taken from Justtanont 

intermediate mass stars for only a relatively short period 
(~ 10 4 yr). From our calculation, the gas-to-dust mass ra- 
tio ranges from 50 to 180, with a mean value of 110. This 
is in agreement with the value of 160 from Knapp ( 1985). 

3.3. Correlation with the optical depths 

From fitting the SED to the observations, we divide our 
sample into two main categories. The first group con- 
tains stars which do not show the ice absorption at 3.1/im 
(Fig.GJ. We modelled these stars using pure silicate grains. 
The second group of stars do show ice absorption from 
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symbols are the same as in Fig. [21 These stars have much 
to the dust grains to fit the SED (long-dashed line). The 
& Tielens lfTM?j) . 

SIRCA and/or CGS4 observations. Their SEDs can be 
described by the two components of dust, i.e., warm in- 
ner silicate grains and composite grains where the tem- 
perature drops below 100 K. This group can be further 
subdivided into those which are well fitted by the same 
silicate grains used for group 1 (Fig. |2J), and those which 
show near-IR deficiency compared to the model (Fig. |3J) . 
For simplicity, we refer to groups 1, 2a and 2b. Due to 
the lack of near-IR spectrum of the stars in group 2a, we 
cannot firmly conclude whether the silicate dust for these 
stars should be different from group 2b. We will discuss 
this point on the following section. 
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Table 3. Input parameters for the modelling and the derived silicate and water-ice optical depths. 
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14.72 
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a Sevenster ((2002 ) 
b David et al. ftMfl 
c Kemper et al. lT?003l 
te Lintel Hekkert (1989) 



We find a strong correlation between the derived mass 
loss rate and the depth of the 9.7/im feature (FigQJ. We 
assume here that the interstellar contribution is negligible 
compared to the circumstellar contribution. The depen- 
dence is not surprising as the depth of the silicate feature 
is a direct measure of the amount of dust along the line 
of sight. The relationship between the measurable depth, 
simply the ratio of the flux at the bottom of the silicate 
feature and the continuum over the feature at the same 
wavelength can be expressed as 

log M = -3.1 - 3.4 (1) 

-^cont 

Here, we defined a continuum flux at 9.7//m as a lin- 
ear interpolation of the fluxes between wavelength of 8 
and 12/im over the silicate absorption feature. The un- 
certainty of the relationship is a factor of 0.2 on the log 
scale (3cr) which translates to a factor of 1.6 in the mass 
loss rate. A similar relationship between mass loss rate 
and the strength of the silicate feature has been shown to 
be the case for optically thin AGB stars and supergiants 
using the silicate emission (Skinner & Whitmore |T998). 
Our derived relation simply extends the estimation of the 
mass loss rate into the very optically thick regime, i.e., 
M > 1(T 4 M Q yr" 1 . 

In Fig. |SJ we show the plot of the derived mass loss 
rate and the derived optical depth at 3.1/im, T3.1 from the 
SED fitting. Note that the value of the optical depth in 



Table |2| is the total optical depth due to both the pure sil- 
icate grains and silicate core-ice mantle grains hence the 
values are relatively high. The correlation is very strong, 
regardless of whether the star shows the ice absorption 
band. We find the same fit even if the stars in group 1 
(those with no water-ice) are not taken into account. This 
confirms the correlation seen earlier by Meyer et al. jl998) 
for a smaller sample of star. They stipulated a direct trend 
that a star with a high mass loss rate shows deeper ice ab- 
sorption and hence higher column density of water-ice. It 
is somewhat interesting that stars in group 1 also show 
the same trend of higher T3.1 with the mass loss rate. This 
reflects the very high silicate dust absorption efficiency 
without invoking other dust type. Note here that the op- 
tical depth in table [3] and Fig. [S] is the total optical depth 
integrated along the line-of-sight and includes the contri- 
bution of both silicate and water-ice. Most of the opacity 
is due to the warm silicate dust which rises with decreas- 
ing wavelength hence giving the correlation seen for group 
1 stars. As can be seen in Fig. not all stars with very 
thick circumstellar shell show evidence of ice formation. 
There is also a slight hint that the relationship may taper 
off towards higher optical depth and mass loss rate. 
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Fig. 4. A plot of the derived mass loss rates against the 
depth of the 9.7/im silicate absorption. The line represents 
the best fit to the data. The filled and opened symbols are 
stars with and without water-ice absorption, respectively. 
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Fig. 6. The ISO-SWS continuum divided spectrum of 
OH26.5+0.6 (solid line) compared with the calculated 
absorption spectrum for the gas temperature of 300K 
(dashed line) and 500K (dotted line). 














































o A 












O 






>> 


















AA 






o 




A 




A 








O 








CD 




A 








at 




A 












A 








w 
w 


o 










Q 




A 




















£/] 


o 










W 


o 






















>: 


o 












o 












A 












o 












o 










io- 4 













2 4 6 8 

T 3.1 



Fig. 5. A plot of the mass loss rates against the optical 
depths at 3.1/zm for stars with (solid triangle) and with 
no (open circle) detected water-ice absorption band. 

3.4. Ice condensation 

It seems that even though the condition for ice formation 
is a highly optically thick envelope (high density) so water- 
ice can be shielded from the central stellar radiation field 
(low temperature), not all the stars in our sample show 



the ice absorption band at 3.1/im. From simply knowing 
the mass loss rate, it is not possible to predict which star 
has the icy mantle grains. However, it is true that for 
those which have water-ice signature, the mass loss rate 
correlates with the optical depth, and hence the ice column 
density, iVj ce , at 3.1/xm via (d'Hendecourt & Allamandola 

N ice = T 3 .iis/A (2) 

where v and A are the full-width-half maximum and the 
absorbance of the water-ice band at 3.1/mi, respectively. 

It is a puzzle why some stars form water-ice and not 
others. One possible explanation is the water vapour abun- 
dance in the photosphere and circumstellar envelope. We 
therefore investigate the existing spectra of these stars. It 
is unfortunate that the near-IR spectra of the sample stars 
are not good enough to see the stretching bands of water 
at 3/im from the available ISO-SWS data as the signal-to- 
noise ratio in these spectra is very low. The bending band 
around 6/xm were then checked to see whether there are 
corresponding lines. In this region, there are expected to 
be strong lines which are sensitive to the temperature of 
the gas, i.e., at 6.59, 6.63 and 6.83/im. The ratio of the first 
two lines are a good temperature indicator (Fig. EJ. These 
lines are seen in the continuum divided spectra of the sam- 
ple. The ratio indicates a temperature in excess of 500K. 
However, the expected strong feature at 6.83/im is not 
distinctly detected, pointing to a lower (T<300K) tem- 
perature gas. Other weaker features cannot be confirmed 
so the detection of gaseous water remains tentative. The 
overlap between our sample and the SWS spectra happens 
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to be those stars which show water-ice signature. The lack 
of the detection of water vapour lines may not be surpris- 
ing taking into account that much of the available water 
freezes out onto grains. In this scenario, gaseous water 
does not provide a significant opacity source in the near- 
IR. 

In order for water-ice to condense, the density has to 
be high and the grains must be cool enough for the water 
to condense onto the existing silicate grains. This requires 
that the dense circumstellar envelope has to be of an order 
of a few 10 16 cm as the dust temperature falls off as r~ a4 , 
according to our radiative transfer calculation. This is in 
accord with the calculation by Dijkstra ct al. (2003a) who 
estimated the radius of the water-ice condensation zone 
started from around (1-2) 10 3 AU, the exact distance de- 
pending on dust grain radius. However, we are not able 
to measure the radius of the superwind and it is possible 
that in some of the stars, the onset of the superwind is 
very recent that the outer edge of the shell has no time 
to cool down below the water-ice condensation tempera- 
ture. In such a case, we do not expect to detect the 3.1/im 
absorption. 

4. Near-IR opacity in icy envelopes 

It is clear that in stars with a clear signature of water-ice 
absorption in the envelope (group 2b) , the near-IR is sup- 
pressed compared to those with a small amount of (group 
2a), or no (group 1), water-ice. The spectrum shortward 
of 5/im cannot be explained by the presence of the same 
type of silicate or the water-ice. In fitting the spectrum of 
these stars, an adhoc near-IR opacity is included in the sil- 
icate grains (see Fig[3J|. As we did not see a clear evidence 
for the gaseous water in the spectra of these objects, we 
search for an explanation in terms of dusty components 
to provide the needed opacity. 

It is known that very optically thick envelopes around 
extreme OH/IR stars exhibit crystalline silicate feature in 
the far-IR (e.g., Cami et al.. 11998). However, the near-IR 
opacity of the crystalline silicate is much lower than that 
of the amorphous grains (Koike et al.. I1993|) . This option 
cannot then be used to explain the near-IR deficiency seen. 

Kemper et al. I|20021 reported on an extreme OH/IR 
star, OH127. 8+0.0, which also has a near-IR deficiency. 
Although this star show an evidence for water-ice in its 
spectrum, simply including it in the modelling could not 
explain the observed deficiency. They, instead, proposed 
metallic iron as being the absorber in order to fit the 
observed spectrum. It remains curious that the lack of 
near-IR flux is seen to be accompanied by the water-ice 
absorption, as is the case in our sample stars. It is possible 
that when the water-ice condensed onto the silicate grains, 
some other materials are incorporated in the matrix lead- 
ing to the extra absorption in the near-IR. However, the 
link between the condensation of water-ice and iron grains 
is unclear, if the extra opacity in the near-IR is due to iron 
dust. Thermodynamically, the two are not linked due to 
the difference in the condensation temperatures. Iron is ex- 



pected condense close to the central star (Gail & Sedlmayr 
1986). By removing iron from the gas phase, it naturally 
leads to the explanation of magnesium-rich silicate grains. 
Many OH/IR stars have been found to have the signature 
of fosterite, a pure Mg-silicate (Waters & Molster ©99 ) . 
Although it is viable that the near-IR opacity is due to 
iron, it remains to be seen why iron condenses out (close 
to the star) at the same time water-ice condenses onto 
cold silicate grains farther out. 

Alternatively this could simply be due to the way these 
objects are viewed. It is possible that these stars have a 
disk which is where the water-ice resides. When this disk is 
along the line of sight (edge-on disk), the absorption is at 
its maximum and hence the the near-IR radiation is more 
efficiently absorbed by the disk while stars which the disk 
is perpendicular to the line-of-sight will show much less 
absorption. Such models have been demonstrated by Ueta 
& Meixner l)2()03|l who calculated a 2-D dust radiative 
transfer model. When viewed edge-on, the near-IR flux 
drops relative to when the disk is viewed face-on. If this 
is the case, it implies that already at this stage of the 
AGB evolution, the star can be seen to depart from having 
spherically symmetric outflow. The presence of a disk will 
enhance further mass loss in an axisymmetric manner. It 
is known that most of the observed post- AGB stars show 
axisymmetric morphology (Ueta et al.GiOOO). 



5. Enhanced mass loss - superwind 

To investigate further, we check the OH 1612 MHz maser 
flux for each object. Here, there is no real correlation be- 
tween the maser luminosity and the ice condensation as 
the stars with and without ice are intermingled. It also 
does not follow the expected trend between the OH maser 
luminosity and the mass loss rate, namely M/v oc Lq^ 
(Baud & Habing 1983). For an expected OH luminosity, 
the derived dust mass loss rate is much higher than ex- 
pected. This deviation can be understood in the context 
of the superwind. These stars have most likely recently 
undergone an increase in the mass loss rate which has not 
have time to propagate out to where the OH maser emis- 
sion originates from. Typically, OH maser scale is of an 
order of 10 17 cm which reflects the mass loss rate of about 
2 000 years ago. The discrepancy between the dynamical 
and the expected mass loss rate derived from the OH lu- 
minosity reflects the time scale of the superwind of less 
than about (1-2) 10 3 years old. An enhanced mass loss 
rate varies by about a factor of 2-6 compared to those 
derived from the OH maser. 

This is further confirmed when looking at available CO 
observations of our sample. Only very few have been ob- 
served and detected (Table ^J. Their CO mass loss rate is 
significantly lower than the derived total dynamical mass 
loss rate by a factor of 10 or more. The estimated CO 
mass loss rate is based on the formalism by Olofsson et 
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Table 4. Calculated gas mass loss rates from dynamical 
calculation (present work) and those derived from OH lu- 
minosity and CO observations. 



Star 


thio "work" 


OH 


CO 


TR AS 1 8059 


1.1E-4 


6 3F-05" 




TO A Q 1 SnQ9 


1 1P A 






TR A S 1 nn 


9 IF A 

Z . 1 Hi- ^± 


'-t.yjlli-yJO 




OH 16.1 


2.9E-4 


4.7E-05 a 




TR AS 1 89^7 


1.5E-4 






OH 21.5 


2.6E-4 


8.1E-05™ 




TRAS 18327 


1.2E-4 


5 0F-05 a 

'J . V/l — 1 Vj!j 






2.6E-4 


q DF D^™ 


8 DF Pi d 


TRAS 18361 


1.7E-4 


5.1E-05" 




OH 30.7 


2.6E-4 


4.5E-05" 




OH 30.1 


2.2E-4 


7.5E-05 a 


1.2E-5 e 


IRAS 18488 


3.6E-4 


6.5E-05" 




OH 32 8 

V / 1 1 tl^j<U 


3.1E-4 


1.6E-04 6 


3.1E-5 e 


TRAS 18588 


1.3E-4 


2 5F-05™ 




OH 39.9 


1.6E-4 


6.6E-05 a 




IRAS 19067 


2.7E-4 


4.9E-05 a 




IRAS 22177 


1.3E-4 


1.0E-04 b 


6.5E-5 d 


AFGL 5379 


2.0E-4 


5.3E-05 a 


3.4E-6 / 


OH 13.1 


2.4E-4 


8.3E-05 6 




IRAS 18310 


2.4E-4 


3.8E-05 C 




OH 51.8 


1.8E-4 


8.8E-05 C 




OH 65.7 


3.2E-4 


2.6E-05 C 




OH 75.3 


3.0E-4 


5.9E-05 C 





a Sevenster (I2002l> 
6 te Lintel Hekkert et al. l|1989|l 
c David et al. (TlMfl l 
d J=3-2 from Kemper et al 
e J=2-l from Heske et al. fiMty 
1 J=2-l from Loup et al. fTMfll 



al. and Gronewegen et al. Ijl999j) who extended the 

equation for CO J=3-2. 



M = 



1.4 T MB v 2 D 2 e 2 
210 19 /cos(J) 



(3) 



where Tmb is the main beam temperature, v is the termi- 
nal velocity in km s , D is the distance in pc, 9 is the 
FWHM of the beam, fco is the CO abundance (taken to 
be 5 10 -4 here) and s(J) is a constant, depending on the 
rotational transition. The CO mass loss rate can be seen to 
be lower than the OH mass loss rate which supports that 
CO transitions probe the much earlier mass loss episode 
than that of OH and the infrared. The disparity between 
the dynamical mass loss rate and those derived from CO 
observations lends a very strong support to the superwind 
idea. 

By comparing the OH mass loss rate and the dynami- 
cal mass loss rate (Table QJ, we can determine the upper 
limit of the lifetime of the superwind of 2 000 yrs. This 
gives enough time for the dust to cool below the water-ice 
condensation temperature. What is not clear is the ac- 
tual time since the superwind started in each object. In 
some stars, the superwind might have just started a few 
hundred years ago hence grains are too hot for water con- 



densation. However, Justtanont et al. (1996) showed that 
OH26.5+0.6 entered the superwind phase some 200 years 
ago and still the spectrum shows the evidence for water- 
ice (Meyer et al 1998 ). This means that the high density of 
the superwind can effectively shield the stellar radiation, 
allowing silicate grains to be cool enough for water-ice 
condensation. In order to determine the extent of a super- 
wind, high resolution observations in the infrared, or CO 
and modelling will be needed (e.g., Fong et al.. I2002jl . It 
may also be possible to probe the superwind using high 
angular resolution OH maser observations as these stars 
are known to have relatively strong maser emission. 

6. Summary 

We search the circumstellar envelope of extreme OH/IR 
stars for the presence of water-ice in an attempt to check 
the reported correlation between the water-ice column 
density and the mass loss rate (Meyer et al I1998|l . We 
found that such a correlation exists in our sample stars 
which display the 3.1/im water-ice absorption. However, 
not all the stars in our study show the water-ice band. 
Furthermore, from the depth of the silicate 9.7/zm band 
which is an indication of the dust mass loss rate, it is not 
possible to predict which of the circumstellar shells should 
display the water-ice band. Stars which show a distinct 
water-ice absorption band also shows a near-IR deficiency 
compared to stars with no water-ice detection. A plausi- 
ble explanation is another dust species which condensed 
at the same time as the water-ice, e.g., iron (Kemper et 
al. [2002 ) . It remains a puzzle how the two species with 
different condensation temperatures are linked. It is also 
possible that the water-ice is present in a disk around the 
central star and when viewed along the line-of-sight, gives 
higher extinction than at an inclined angle. 

We found a simple estimate of the mass loss rate from 
the ratio of the depth of the silicate absorption feature at 
9.7/im and the continuum flux. The main uncertainty in 
our estimated dust mass loss rate is the distance to the 
stars which in many of the cases here is unknown. 

Comparison of our derived dust mass loss rate with the 
expected gas mass loss rate from the OH luminosity and 
CO rotational observations indicates that these extreme 
OH/IR stars have recently undergone an increase in the 
mass loss rate by a factor of > 10. The maximum timescale 
for the superwind is less than the time it takes for the 
mass loss to reach the H 2 photodissociation radius, i.e., 
~ 2 000 yrs. In some cases, this can be much smaller as the 
dynamical mass loss rate is higher than the mass loss rate 
derived from the OH observations. High angular resolution 
observations (< 0.1") will be able to confirm the extent of 
the superwind. 
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